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G\ ■ ABSTRACT 

1^ • We present arcsecond-scale mid-ir photometry (in the 10.5 fim N band and 

I . at 24.8 /im), and low resolution spectra in the N band {R ~ 100) of a candidate 

2 ■ high mass protostellar object (HMPO) in IRAS 18151-1208 and of two HMPO 

^ '. candidates in IRAS 20343+4129, IRS 1 and IRS 3. In addition we present high 
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resolution mid-ir spectra {R ~ 80000) of the two HMPO candidates in IRAS 
20343+4129. These data are fitted with simple models to estimate the masses of 
gas and dust associated with the mid-ir emitting clumps, the column densities of 
overlying absorbing dust and gas, the luminosities of the HMPO candidates, and 
the likely spectral type of the HMPO candidate for which [Ne II] 12.8 /im emission 
was detected (IRAS 20343+4129 IRS 3). We suggest that IRAS 18151-1208 is 
a pre-ultracompact HII region HMPO, IRAS 20343+4129 IRS 1 is an embedded 
young stellar object with the luminosity of a B3 star, and IRAS 20343+4129 IRS 
3 is a B2 ZAMS star that has formed an ultracompact HII region and disrupted 
its natal envelope. 

Subject headings: stars: formation — stars: pre-main sequence — circumstellar 
matter — ISM: jets and outflows — HII regions — infrared: ISM 



Introduction 



Many open questions in high-mass star formation are related to the evolution of cir- 
cumstellar envelopes, accretion disks, and jets from high-mass protostellar objects (HM- 
POs). HMPOs are often bright sources in mid-ir continuum, but only in a few recent cases 
have images suggested specific structures such as disks, jets, or warm o utflow cavity walls 
dSridharan et ahlboosl : IPe Buizer fc Minier|[2005l : IPe Buizer Ibood . l2007ah . Mid-ir ionic lines 
like [Ne II] and [S IV] have been used to map compact and ultracompact HII r egions (UC 



HII) and photodissociation regions, and to study their structu r e and excitation (ILacy et al. 



19821 : lOkamoto et al.ll200ll : IZhu et al.ll2005l : iKassis et al.ll2002l . 120061 1 . but there is a lack of 
observations of HMPOs. A remarkable feature of high-mass star formation is that HM- 
POs, defined as actively accreting mass, can begin nuclear fusion and hence also be rapidly 
evolving massive young stella r objects (MYSOs) that have already formed hypercompact 
or ultracompact HII regions (iBeuther et al.ll2006l ). This feature raises a possibility of de- 
termining the spectral type of an MYSO through the ionic lines excitation, or from the 
number of ionizing photons required for its observed centimeter continuum flux, separately 
from estimating its luminosity and spectral type from infrared emission. However there is 
also the possibility that the ionization is collisionally excited by a jet. It may be possible to 
distinguish between the two cases, depending on the Doppler velocities, the morphology of 
ionized gas, and the ratio of the required flux of ionizing photons to total luminosity. 

Hoping to enlarge the sample of HMPOs that could be studied in detail despite po- 
tential limitations, in 2003 we made mid-ir observations on the IRTF of about a third of 



the survey of 69 HMPO candidates presented by ISridharan et al.l ( 120021 ) and iBeuther et al. 
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( l2002al ). We chose objects from their survey that appeared compact and/or bright in the 
MSX su rvey, and found that about 80% of them were unresolved or marginally resolved by 
MIRSI (jPeutsch et al.l |2002| ) on the IRTF in the broad N band at 10.5 fim and in a nar- 
row band filter at 24.8 /im. In addition, we obtained MIRSI grism low- resolution spectra 
(R ^ 100) of teii of them in the N band. In 2006 on Gemini North, we obtained TEXES 
(ILacy et al.ll2002l ) high- resolution spectra (R ~ 80000) of two HMPO candidates for which 
we had grism spectra. In this paper we present spectra and photometry of three candidate 
HMPOs including the two with TEXES spectra: IRAS 18151-1208, IRAS 20343+4129 IRS 
1, and IRAS 20343+4129 IRS 3 (hereafter, 18151, 20343 IRS 1, and 20343 IRS 3). We will 
demonstrate that mid-ir emission from the dust and gas near the HMPO candidates (where 
it is strongly heated) can be used as a useful probe of temperatures, masses, and luminosi- 
ties, using simple isothermal clump models, even if each component (envelope, disk, jet, or 
cavity wall) is not resolved. In combination with observations cited below, we are able to 
use our new data to infer the nature of each candidate candidate HMPO (e.g. pre-UCHII 
region HMPO, ZAMS B2 star). 

The objects c hosen are near the cen ters of complex, large-scale massive molecular out- 
flows mapped by iBeuther et al.l (l2002bl ). 20343 has an apparent larg e-scale N-S outflow 



whose red and blue lobes are both extended E-W (IBeuther et al.ll2002bl ). but IRS 1 also has 



fPalau et al. 


2007) 


Kumar et al. 


2002) 



20021 ). Two of them (1815 1 and 20343 IRS 3) were observed to hav e 0.5 and 



1.8 mJy 3.6 cm emission, respectively (iCarral et al.l 119991 : iSridharan et al.ll2002l ). We ob- 
served 20343 IRS 1 and IRS 3 with TEXES on Gemini North based on the 3.6 cm and H2 
observations, with the goal of studying the role of ionized gas in them. 

The 10 /im grism spectral shapes of the HMPO candidates fall into three classes: those 
with deep silicate absorption; those with moderate sihcate absorption and an apparent peak 
at about 8.5 fim; and those without an apparent silicate absorption feature but with contin- 
uum rising monotonically from short to long wavelengths (Campbell et al. 20 07, in prepara- 



tion) . Examples of these shapes can be seen in the UC HII spectra presented by lFaison et al 
( 119981 ). Since the HMPO candidates were chosen based on IRAS colors similar to UC HII 



regions (ISridharan et al.ll2002l ). one would expect the HMPO candidates to have similar 10 
/im spectra. IR AS 20343+4129, was observed with the IRAS LRS and has a silicate ab- 
sorption feature (jVolk et al.lll99ll ). The three objects presented here include an example of 
each of the three classes of grism spectral shapes. Two of the ten c andidate HMPO grism 
spectra show strong [Ne II] lines, IRAS 18247-1147 and 18530+0215. ISridharan et"aD J2OO2I ) 
reported relatively strong 3.6 cm fluxes of 47 and 311 mJy, respectively, for them. We did 
not observe them with TEXES on Gemini in order to focus on the earliest possible HMPO 
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stage associated with ionized gas. 



Deriving physical information from the continuum spectra is difficult because the actual 
geometry of the dust distribution is unknown, except that the sizes of the N band and 24.8 /xm 
images limit the extent of the emitting structures. Experience with one -dimensional radiative 



trans f er mo dels of spectral energy distributions from UC HII regions (jCampbell et al.lll995 



2OOOI . I2OO4J ). and inspection of the new spectra themselves suggest that there are ranges 
of temper atures in the ern i tting re gions . However, the two- di me nsion al radiative transfer 
models of IPe Buizer eTaP J2005bh and Iwhitnev et aP J2003al lbl. l2004h show that orienta- 
tion of flattened envelopes and outflow cavities dramatically affects the depth of the silicate 
absorption feature, as does emission and a bsorption by individual c lumps in a clumpy dust 
cloud in the three-dimensional models of Indebetouw et al. J2006h . The recent observa- 
tions of extended and complex near- and mid -ir emission from HMPOs (jPe Buizer &: Minier 
2OO5I : iDe Buizer II2006I : ISridharan et al.ll2005l ) also indicate that one-dimensional models are 
unrealistic. Nevertheless, a simple three part model will allow us to derive approximate tem- 
peratures, column densities, and masses of different dust components. The first component 
is hot dust that could be (part of) a relatively compact disk, or a clump very near the HMPO 
candidate; the second is warm dust that could be a more extended (part of a) disk, a clump 
of dust further out, or perhaps the inner wall of of an outfiow cavity; and the third is cold 
dust in an outer envelope that creates the silicate absorption feature. 

Deriving information from the ionic lines of an UC HII region is in principle straight- 
forward. The line fiuxes can be corrected for local extinction based on the continuum models 
discussed above, and then the ratio of [Ne Ilj to [S IV] fiuxes can b e used to determine the ex- 
citing star's temperature (ILacy et al.lll982l : lOkamoto et al.ll200ll ). The numerical simulation 



code cloudy (http://www.nublado.org/) can also be used to deduce the star's temperature 
and luminosity by matching the simulated intensity and spatial extent of free-free, [Ne II], 
and [S IV] emission to the observations. The star's temperature can then be compared to 
that of the spectral type deduced from the cm continuum fiux. In addition, comparison of 
Doppler velocities of the ionic lines to those of molecular lines can be used to indicate if the 
gas is in an UC HII region or a jet. 



Observations and Reduction 



2.1. MIRSI Images on the IRTF 



MIRSI, the Mid-InfraRed Spectrometer and Imager (IDeutsch et al.ll2002l ). was operated 
remotely from Colby College and The Center for Astrophysics on the IRTF for these ob- 
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servations. MIRSI has an imaging field of view of 85 x 65" on the IRTF, with diffraction 
limited performance (0.8" at 10/im), and plate scales of 0"2689/pixel in right ascension and 
O'.'2635/pixel in declination. The objects presented here were imaged on 2003 September 13. 
From a trial observation of HMPO candidates during a 2002 engineering run, we expected 
the sources to be unresolved or marginally resolved. Conventional chopping and nodding 
was used with chopper throw 25" NS and telescope nod 25" EW, so that all chop-nod images 
were on the array. Five dithered integrations with total on-source time of 240 seconds were 
taken through a 5 fim wide N band filter centered at 10.5 fim and then through a 7.9% wide 
filter centered at 24.8 /im of the 18151 and 20343 fields. Both 20343 IRSl and IRS3 were 
present in the latter field. The dither pattern was a central position with ±5" offsets in right 
ascension and declination. 



^Excess non-statistical electronic noise was removed with a custom IDL procedure (iKassis et al 

2OO3)0 Chop-nod addition and subtraction and fiat fielding were done in the usual way. Neg- 
ative chop-nod images were inverted, and all individual images from the dithered chop-nod 
sets were combined in IRAF. 

Positions were determined relative to those of calibration stars without special care for 
precision astrometry, so the positions are limited by the inherent offsetting accuracy and 
stability of the IRTF on that night to several arcseconds. Since 20343 IRS 1 and IRS 3 
were observed simultaneously in each individual frame, their relative positions should be 
accurate to better than two pixels (0"5). The relative positions of the r nid-ir sources at the 
two wavelengths agree to 0"1 and they agree with the K band sources of iKumar et al.l (|2002| ) 
to 0'.'5, although it is not known if the mid-ir centers fall exactly on t he K band centers . 
We assume that 18151 is coincident with the peak of 3.6 cm emission (jCarral et al.lll999l ). 
Positions are shown in Table 1 for the sources. 

Simple photometry was performed using the IRAF task imexamine. The K stars 7 
Aql and 7 Dra were used for point spread functions (PSFs) and photometric flux density 
calibrations. N band magnitudes for calibration were taken from the list of bright infrared 
standard stars on the IRTF web site. We used the N band effe c tive w avelength (10.47 fim) 
and the flux density for zero magnitude as given by iTokunagal (I2OOOI ). For the MI RSI 24.8 
/ im fil ter, we used a zero magnitude flux density calculated by the formula given bv lEngelke 



19921) shifted to agree with the N band zero magnitude flux density given by iTokunaga 
(I2OOOI ). We assumed that the 24.8 /im magnitude would equal the magnitude in the 20.13 fim 



^The MIRSI camera is divided into sixteen, 20-colunin channels, with each column containing 240 pixels. 
For each image frame, the pixels in all channels displayed a distinctive noise pattern common to all that was 
determined and removed. Each channel had a different median offset that was also removed. 
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Q band (ICohen. et al.lll995l ). This procedure gave agreement with the IRAS color corrected 



25 fim flux density of 27.0 Jy for 7 Dra within 1%. We did not detern aine atmospheric 



extinction coefficients, and used the mean values for N and Q given by iKrisciunas et al. 



( Il987l ). N band flux densities appear to be reproducible to ±10% and to have a similar 
systematic accuracy, but the 24.8 fim reproducibility and systematic accuracy is probably 
±30% due to the variations in the atmospheric extinction over the course of the night. The 
flux densities are shown in Table 1, including the values used for 7 Aql. 

Images of 18151 and 20343 IRS 1 appear to be unresolved at both wavelengths with in- 
dications of the first diffraction ring. 20343 IRS 3 is marginally resolved at both wavelengths 
with FWHM a few tens of per cent larger than the PSF. In addition, there is extended diffuse 
emission surrounding the central peak, especially in N band (Fig. la). The extended nature 
of the central peak of IRS 3 is clearer at 24.8 /im, and can be seen in Fig. lb. FWHM values 
of Gaussian fits to the profiles by IRAF are shown in Table 1. Our survey observations 
and "pipeline" image processing did not emphasize the highest possible S/N for source and 
PSF profiles intended for image deconvolution. We characterize the PSFs by values shown 
in Table 1 that are averages of Gaussian FWHM values and their standard deviations from 
calibration star observations over the course of the night. The standard deviation in the PSF 
could be due to small errors in shifting individual images before averaging and/or varying 
seeing in the unstable atmosphere. 



2.2. MIRSI Grism Spectra on the IRTF 

Grism spectra of the sources were obtained on 2003 September 14. The MIRSI grism 
was used with an N band prefilter, cutting off the spectra just longward of [Ne II] 12.8 fim. 
The slit was oriented NS, and the system can be thought of as a long slit spectrometer. On 
E. ToUestrup's suggestion, we arranged both chop and nod to be NS so that four spectra 
were recorded in each nominal 30 second integration camera frame, giving total on-source 
time of 120 seconds per camera frameEl We experienced considerable difficulty centering the 
source on the slit, and only one of the chop- nod images could be well centered on the slit at 
best. For each source, multiple nominal 30 second integration frames were taken, although 
not all could be used due to the source drifting off the slit or excessive noise. 



^MIRSI's electronics adds multiple chip frames, in this case of 700 ms duration, recorded as extensions 
in fits files. Separate fits extensions are recorded for each chop and nod position, and are then appropriately 
added or subtracted to make up a camera frame. The nominal integration time would be the on-source time 
if the chopping and nodding were off the chip. 
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The initial spectral processing involved the following steps: (1) For each chop-nod sub- 
tracted camera frame, the electronic pattern noise was determined from the channel of the 
chip blanked by the N band filter's long wavelength cutoff, and subtracted. (2) The frames 
were flat-field corrected using a dome fiat from which a dark frame had been subtracted. (3) 
All camera frames were averaged into a single combined frame with the IRAF task combine. 
(4) Copies of the combined frame were inverted and shifted as necessary, and then all four 
chop-nod spectra were averaged to a single combined spectrum using the combine task. (5) 
A bad pixel mask was applied to the single combined spectrum. (6) Each channel's dark 
sky baseline surrounding the combined spectrum was inspected for an offset. If an offset 
was found, the median level was shifted to obtain zero baseline outside the spectrum. This 
processing proved capable of recovering spectra where there was no signal visible in raw data 
frames. 18151 had strong signal visible, but both sources in 20343 appeared only barely vis- 
ible after the first stage of pattern noise correction. IRSl suffered from a noisy portion in its 
two best frames between 10 and 11 /zm. 22% of the pixels in the spectrum in this wavelength 
band had excessive noise and were zeroed. The final flux densities were corrected for the 
zeroed pixels. 

For the sources here, fj, Cep, an M supergiant, was used to calibrate the spectra. M 
supergiants are variable, and less desirable than the K giants used as photometric standards, 
but n Cep gave an extremely high S/N calibration spectrum. The fi Cep grism spectra were 
processed in the same way as the candidate HMPO spectra. The IRAF task apall was used 
to extract the candidate HMPO spectra and the fi Cep spectrum. A linear grism dispersion 
relation was created using the telluric lines near 10 fim and the [Ne II] 12.8 fim line. The 
atmospheric transmission spectrum was taken from the UKIRT web site, convolved with the 
grism response, and cross-correlated to the uncahbrated MIRSI fi Cep spectrum to determine 
one wavelength-pixel point. The [Ne II] 12.8 fim line was used from the IRAS 18247-1147 
spectum in which it was detected for the second point. 

In order to use our fi Cep spectrum for flux calibration, we obtained ISO SWS spectra 



from the University o: 



Cal 



gary web sitelf] This site uses a program to extract spectra as 



given by ISloan et al.l (120031 ). Three spectra are available for fi Cep, and they show some 



minor differences. We chose TDT 05602852 for /i Cep becaus e it appeared b est based on 



inspection of separate up and down scans kindly provided by iKraemerl (120061 ). A custom 
IDL procedure was used to calibrate the candidate HMPO spectra to the ISO /i Cep flux 
densities. It multiplied each uncalibrated candidate HMPO spectral point by the ratio of 
the the ISO calibrated spectral point (convolved to the same resolution as the MIRSI grism) 



3 



http: / / www.iras.ucalgary.ca/satellites/~volk/getswsspec_plot.html 
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divided by the uncalibrated MIRSI point for jj Cep. This process gave us a prehminary 
cahbrated spectrum. However, we did not know how well the sources were centered on the 
slit, so a final step was made, following a suggestion by E. ToUestrup. For each preliminary 
calibrated spectrum, the flux densities. Fx, were summed over the spectrum, and the sum 
compared to the flux, F, as measured photometrically through the N band fllter. Their 
ratio was used as a correction factor for preliminary calibrated spectrum. Flux densities in 
Janskys, F^, were calculated from Fx for presentation in this paper. The calibrated spectra 
are shown in Figure 2. The 24.8 /im flux densities in Table 1 compared to the spectra indicate 
that the SEDs all rise signiflcantly with increasing wavelength. 



2.3. TEXES Spectra on the Gemini North 

High spectral resolution observations of 20343 IRS 1 and IRS 3 were made during the 
Texas Echelon Cross Echelle (TEXES) Demonstration Science run on the Gemini North 8-m 
telescope in July 2006 as part of the program GN-2006A-DS-2. TEXES is a cross-dispersed 
mid-infrared (5-2 5 /xm) spectrograph capable of 0'.'4 and 3 km s~^ resolution on Gemini 



( iLacy et al.l |2002| ) . All data from the Demonstration Science run are publicly available at 



http: / / archive.gemini.edu 



The TEXES candidate HMPO observations were made in the TEXES hi-med spectro- 
scopic mode with a 0'.'5 slit giving ~ 4 km s~^ resolution. The slit length was ~ 4", oriented 
EW, with O'.'15/pixel sampling along the slit. Two observing modes were used: nod mode, 
in which the telescope was nodded at ~ 0.1 Hz to move the source by I'.'T along the TEXES 
entrance slit, and scan mode, in which the telescope was moved south across the sky in 0'.'25 
steps without nodding to map the object. 

Two spectral regions were observed, centered at [Ne II] (12.8 fim) and at [S IV] (10.5 
fim). The spectral coverage at each setting was ~ 0.6%. Spectral calibration was determined 
from sky emission lines, and is accurate to 1 km s"^. Intensity calibration was obtained 
from observations of an ambient temperature blackbody, and for scan mode observations is 
accurate to ~ 20%. In nod mode, there is an additional uncertainty due to the unknown 
fraction of the flux outside of the 0'.'5 slit. For compact sources this is a random uncertainty 
due to seeing and guiding, whereas for extended sources it causes a systematic underestimate 
of the flux. In addition, for both nodded and mapping observations the sky subtraction 
procedure introduces an uncertainty due to the possibility of emission in the sky position. 

20343 IRS 1 was observed only in the nod mode at the [Ne II] setting. Its 12.8 /im 
continuum was clearly detected, with 11.5 Jy measured through the 0'.'5 slit. This flux is in 
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agreement with the MIRSI grism spectrum. No scan map was made, but the source's extent 
along the sht was less than 0'.'5. There was no evidence of the [Nc II] line in the spectrum, 
which is shown in Figure 3. The equivalent width for a narrow (< 10 km s~^) line was less 
than 2.0 x 10^^ cm~^, indicating that the line fiux from a O'.'S region was less than 6 x 10^^^ 
W m~^ {2a uncertainties). We note that the line fiux could be greater than our limit if the 
line source is more extended than the continuum source or if the line is broad. A 100 km 
wide line would have to be 2 x 10~^^ W m~^ to be detected. In addition, the line would 
have been missed entirely if it were at Vlsr — -20 to -50 km s~^, which fell between TEXES 
grating orders, but this would require a velocity shift of >30 km from the molecular 
cloud velocity of +11 km s~^. 

20343 IRS 3 was observed at the [Ne II] setting in both nod mode and scan mode. The 
nod mode observations have higher signal to noise ratio, but have increased flux calibration 
uncertainty. The nod mode spectrum and the scan-mode continuum and line maps are shown 
in Figures 3, 4, and 5. The 12.8 //m continuum flux, derived from a sum over a 1'.'35 x 1'.'25 
region where flux is apparent in the map, is 2.3 Jy. This is about 40% of the flux in the 
MIRSI grism spectrum, suggesting that some extended emission was missed. For the nod 
mode observations, the nod throw was I'.'T EW, and for the scans the sky background were 
taken from positions I'.'S north and south of the peak. Consequently, emission on a scale 
> l'/5 would have been missed. The continuum source appears extended NS in the TEXES 
map, and possibly double-peaked. Its extent is roughly consistent with the size of the MIRSI 
10.5 iim N band image (Table 1; §3.1). The more extended 24.8 iim emission is elongated 
along a NW-SE axis rather than NS (Figure 1; Table 1; §3.1). [Ne II] was clearly detected 
with a spectrally resolved line-to-continuum ratio of ~ 10. The [Nc II] distribution appeared 
more point-like than the continuum distribution, and is located on the northern end of the 
extended continuum source (see Figures 4 and 5). The line is centered at Vlsr = 15.7 ± 1 
km s~^ with an observed FWHM ~ 8 km s~^. The line flux from the map is 1.0 x 10~^^ W 
m~^. 

We attempted to observe IRS 3 in the [S IV] setting, but failed to detect either continuum 
or line emission. The MIRSI grism spectrum indicates that the 10.5 fim continuum is a factor 
of about 2 weaker than the 12.8 /xm continuum, so it should have been detected, although 
not easily since the TEXES sensitivity is a factor of about 2 poorer at 10.5 iim. Given the 
possiblity of a pointing error, especially since this was the first science run for TEXES on 
Gemini, we choose not to quote an upper limit from these observations. 



- 10 - 



3. Models of the Mid-IR Emission 

3.1. Overview of a Simple Dust Continuum Model 

The geometry of the mid-ir emitting dusty clouds is unknown except for constraints 
on their projected diameters. We can make simple three-component models to match our 
spectra and photometric data, and the models will give estimates of the temperatures of 
the dust clumps, the column densities through them, the clump masses, and the mid-ir 
luminosities. Our models should not match the SEDs outside the mid-ir: we would expect 
them to underestimate the observed SED at both ends of the spectrum. 

The model components are (1) a hot component whose size is constrained by the N 
band image, (2) a warm component whose size is constrained by the 24.8 fim image, and 
(3) a cold, pure extinction component that is responsible for the silicate feature. In order to 
understand clearly the relationships of the various model inputs, we review the emission of 
an isothermal, constant density, dusty clump observed through a colder cloud that creates 
extinction but no emission. The observed flux density. Fx, is given by 

Fa = QBx{T){l - e-"=W)e-"'"(^) (1) 

where Q is the solid angle of the clump, Bx(T) is the blackbody intensity, Te(A) is the optical 
depth of the emitting clump, and Ta{\) is the extinction optical depth of the absorbing (and 
scattering) overlying cloud. In general, emissivity, e(A) = (1 — e"'^'^''^^), but for an optically 
thin clump e(A) = Te(A). The optical depth in emission is given by 

Te{X) = K,{X) j Pddl (2) 

Clump 

where -ft'e(A) is the absorption cross section per mass of the emitting dust in cm^ g~^, p^ is 
the mass density of the dust, and the integral is through the emitting clump. The integral 
through the clump can be related to the column density of H nucleons, Nh through 

Pddl = ^pmnNH (3) 

Mg 

Clump 

where ^ is the dust to gas mass ratio, p is the mean molecular weight for the assumed 
abundance ratios, assuming neutral atomic gas instead of H2 in order to follow the convention 



of iDraind (l2003al ) that uses the column density of H nucleons, and mn is the mass of a 



hydrogen atom. Thus the emission optical depth can be related to the column density: 

re(A) = K,i\)^pmHNH (4) 
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We can now clearly see from equation (1) that the observed flux density, Fx is related to 
both the solid angle and the column density, and that the derived column density is strongly 
dependent on the assumed solid angle. For the simple geometry of a constant density, 
end-on cylindrical cloud whose solid angle is estimated from the observed angle, Ofwhmi 
il = jOp^rffM^ the optically thin case 

= lel^HMBx{T)K{X)^^^imHNHe-^'^W (5) 

The assumed projected diameter of the clump is thus clearly a key factor in deriving a rea- 
sonable column density. For optically thick clouds, the exponential function in the expression 
for emissivity makes the column density estimate very strongly dependent on the assumed 
diameter. 

A simple approximate way of estimating source diameters from the IRTF images for 
marginally resolved sources like 20343 IRS 3 (see Table 1) is based on assuming that the 
profiles of the true source, the observed data, and the PSF are all Gaussians. With this 
assumption (that is certainly not correct for an Airy disk), a formal deconvolution would 
give 

= \/^f^P (6) 
where is the true source FWHM, 6^ is the observed data FWHM, and 9p is the PSF 
FWHM. The resultant 9s are 0'.'75 and 1'.'57 in the N and 24.8 fim filters, respectively. In 
this case however, the source extent at the shorter wavelength is more accurately determined 
by TEXES at 12.8 fim on Gemini. Calculated in the same way, the extent is C/.'SG x I'.'l EW 
X NS, with area equivalent to a circular source of 0'.'80. 

For the unresolved sources 18151 and 20343 IRS 1 we can make estimates of the source 
sizes from the IRTF observations based on Gaussian deconvolutions as follows. Let us 
represent 6d in terms of Bp and n standard deviations of the PSF FWHM, ap-. 6d = Op + nap. 
Substituting this expression into Eq(6) gives 

= ^2napep + v?al (7) 

An observation of a source with a true value of Og in which the statistical variation resulted in 
a —nap deviation from Od = 6p + nap would appear to be unresolved with observed 6d = Op- 
Thus we can estimate reasonably likely values for source size 9s from a la deviation by 
substituting n — 1 into Eq (7), and a 3(7 upper limit to a source size by substituting n = 3. 
The resulting values of 9s are larger than ap or 3(7p because of the deconvolution. Table 2 
presents diameters for the hot components and the warm components calculated from the 
Gaussian FWHM values in Table 1 in N band and at 24.8 //m, respectively. They have been 
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converted to AU at the distances of the HMPO candidates for the table. The diameters 
in Table 1 are consistent with those of disk candidates associated with HMPOs cited by 



Cesaroni et al.l (120071 ) . most of which have diameters from 1000-3000 AU. 



The derived column densities depend on th e dust model that specifies K{\). A very 
attractive model is that of lOssenkopf fc Hennind (119941 ) for protostellar cores often referred 
to as 0H5. This dust represents coagulated grains with ice mantles that would be ex- 
pected to form from dust originally in the diffuse interstellar medium during the process of 
molec ular cloud formation. It has b een used to fit submm SEDs of high-mass protostellar 
cores (Ivan der Tak et al.lll999l . l2000l ) and far-ir observations of the UC HII region G34.3+0.2 
( ICampbell et al.ll2004j ). However, when we examined its K{X) behavior around 10 /xm, we 
found that its silicate feature is shifted longward from 9.7 fim and broadened so that it 
does not appear to be compatible with our grism spectra. Models by iDraine J2003al lbl) and 
colleagues fit the shape of our observed sihcate absorption fe ature well. We have chosen 
their 2003 Ry = 5.5 synthetic extinction curve (jDraind l2003al ) for dust in dense clouds in 
the Milky Way. It is accessible on the web and well documented. We have also fit the data 
successfully with Ry = 3.1 dust for diffuse clouds, but Ry = 5.5 dust is appropriate for 
dense clouds and the deduced column densities should be more realistic. The dust model 

5.5 dust to fit 1.2 - 8.0 



has 



Ma 

Ma 



105. iRoman-Zunga et al.l (120071 ) recently found the R 



fim data from the dense core Barnard 59. In retrospect, we found high enough temperatures 
that ice mantles should have evaporated so that the 0H5 dust would not be expected to fit 
our data. 

Fitting the shape of the grism spectrum of 18151 (Figure 2) and the large flux density 
at 24.8 /xm (Table 1) requires a minimum of three components for our models: hot dust 
responsible for the 8 /xm end of the grism spectrum, warm dust for the 13 fim end of 
the spectrum and the 24.8 fim photometry, and cold dust for the depth of the silicate 
absorption. The model calculation has the following input parameters: the temperature of 
the hot component, T^, the diameter of the hot component 9^, the optical depth in emission 
of the hot component at 9.70 fim, Th, the analogous parameters for the warm component, 
Tuj, 9w, and t^, and the extinction optical depth, r^, of the overlying cloud that is too cold 
to emit in the mid-ir. In the mid-ir, the extinction is virtually pure absorption. The model 
uses a spectrum that is the sum of the contributions of both hot and warm components 
each calculated according to Eq(l) to fit the observations. For given optical depths, column 
densities, Nh, can be derived from Eq(4) for each component, using parameters found in 
Draine's web site for the dust model. For end-on cylindrical geometry, the masses of the hot 
and warm components can be derived from A'^^^ for each. Visual extinctions are calculated 
from the extinction cross section per H nucleon at the center of the V band as given in the 
web site and Nh for the cold absorbing component. The SED due to the two dust clumps and 



overlying extinction is calculated to verify that it does not create excessive emission outside 
the mid-ir, and to derive the luminosity of the mid-ir emitting clumps for comparison to 
measurements and estimates of the overall candidate HMPO luminosity 



3.2. Fitting the Continuum Data 

Our goal is to apply our simple three-component model to the grism and 24.8 /im 
photometric data to derive estimates of the temperatures, the column densities, the masses, 
and the luminosities of the emitting dust components, and the column density of the cold 
extinction component. We have no measurements of the angular extents of the components 
for 18151 or 20343 IRS 1 because they were unresolved, although we do have them for 20343 
IRS 3. If we assume incorrectly small sizes for the first two, wc will make over-estimates of 
the column densities. We have chosen to use large size estimates that are consistent with the 
data for these two so that our column densities can be thought of as lower limits to the true 
values. The parameters for our models are shown in Table 3. Table 3 shows sizes based on 
1 a deviation estimates for unresolved sources in Table 2. For 20343 IRS 3, it uses the size 
from the TEXES continuum map for the hot component, and the size based on Gaussian 
deconvolution of the MIRSI 24.8 iim image for the warm component. 

The models were interactively fit to the data. In order to quantify the quality of the 
fit and to aid in choice of extreme values of the parameters consistent with the systematic 
accuracy of the data, we defined a modified statistic. The grism spectra were broken 
into four photometric bands: 8.0 — 9.0/xm, 9.0 — 10.5/im, 10.5 — 12.0/xm, and 12.0 — 13.0;um. 
Fluxes were summed in each band. The 24.8 fim photometry served as a fifth band. We 
defined the modified statistic for each band, i, as 

,Model -^ijData) /q\ 
Xmi ~ p 2 \ ) 

^ i,Data 

where Fj is the flux in band i. The sum of the five terms formed for the model. This 
statistic places the quality of the fit as a fraction of the data value in each band on equal 
footing. Even though it assumes that each of the five bands has equal S/N and has no 
specific statistical interpretation, it is useful for fitting models to the data. For optimizing 
models, wc sought to minimize interactively, initially using graphs of the models plotted 
over the data as guides, rather than using an exhaustive, automated search through the 
parameter space. Such an automated search is not justified by the quality of the data and 
the simplicity of the model. Final models in Table 3 were optimized using values. For 
choosing the upper or lower value of a specific parameter (e.g. Ty^) that might be reasonable, 
we started with the value in an optimum model, and varied that parameter (only) until the 
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band most affected by the parameter (in this case i = 5, 24.8/im) indicated a 30% difference 
(our systematic photometric accuracy), or Xmi — 0.09. On a graph of the data hke Figure 
6, the upper value would cause the model's flux density to fall on the end of the upper 
error bar at 24.8 fim. It was not feasible to vary more than one parameter at a time because 
of the size of the parameter space and the fact that all bands of our data were somewhat 
affected by all of the parameters. The procedure is arbitrary. Larger extreme values for 
a parameter of interest could be found if all others were also varied compensating for the 
changed parameter of interest (e.g. if were the parameter of interest, a larger extreme 
value would result if Th,Tw,Tw, and were varied in addition to Th). 

The best-fit temperatures of the hot and warm components are in the ranges 420-1000 
K and 110-182 K, respectively. It is useful to examine SEDs over a wide range of wavelengths 
in order to understand how the components combine to fit the data. We show the SEDs 
of individual components in Figure 6 for the model of 18151 with a 455K hot component 
whose parameters are given in Table 3. Figures 7 and 8 show models for 20343 IRS 1 and 
IRS 3, and their model parameters are also shown in Table 3. In addition. Figures 6-8 show 
the combined SEDs both with and without the cold components' extinction, and the IRTF 
data. In Figure 6, the 455K hot component shows a strong 9.7 fim silicate emission feature, 
as expected. However the 136K warm component does not show a strong emission feature 
because its Planck spectrum is rapidly rising to long wavelengths and the silicate feature is 
smoothed out because it is optically thick in its center. Absorption features at 9.7 and 18 
/xm due to the cold layer are clear in the total spectrum, and the latter causes a nearly fiat 
portion between 15 and 20 fim. This shape in the mid-ir is not necessarily an artifact of the 
simplicity of a three component model. A similarly shaped SED is shown from the library of 
Monte Carlo radiative transfer models at jhttp: / / caravan.astro.wisc.edu/ protostars/ for the 



low mass protostar IRAS 04368+2557 bv iRobitaille et al.l (120071 ). Our model SED for 18151 



in Figure 6 also clearly shows that these components emit very weakly at both near-ir and 
submm wavelengths. For 18151, a large range in will fit the data apparently because the 
grism spectral range could lie in the Raleigh- Jeans part of the hot component's spectrum 
for high values of Th- The Th = 455 K model is in the lower end of th e range. A model wit h 



Tfi = 1600 K, a commonly assumed dust sublimation temperature (IWhitnev et al 



will f it our data, but would have a large excess over the observed flux at 2.1 fim (IDavis et al. 



2004) 



20041 ). A model with Th = 1000 K is presented in Table 3 that fits both our data and the 



flux at 2.1 fim. 
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3.3. General Results and Discussion for the Continuum Models 

Parameters for the models are presented in Table 3. Two models are given for 18151 
at extreme ends of the range of Th that fit the data well, and the best-fit model is given for 
each of 20343 IRS 1 and IRS 3. The summed values of Xm given in the last column]^ In 
addition, upper and lower extreme values that are consistent with the data to within 30% 
as discussed above are shown. 

There are a number of important aspects of the values derived from the models. The first 
is that the source that appears most deeply embedded, 18151, could have a hot component 
with a sufficiently high temperature so that its mid-ir would be on the Raleigh- Jeans end 
of the spectrum, so that the temperature cannot be determined without measurements in 
shorter bands like K, L or M. However, it is entirely likely that some or all of the flux at 
these shorter bands would come from an additional, hotter dust component than the ones 
whose mid-ir we seek to model. For the other two objects. Table 3 shows moderate ranges 
in temperatures from lower to upper values that are compatible with the best-fit emission 
optical depths. The range of temperatures that could be fit by varying the emission optical 
depths at the same time with only modest increases in would be larger than shown. 
Nevertheless, we feel the temperatures shown are realistic estimates. The ranges in Th and 
are well separated in all of the sources. If each source's emission is from a single clump or 
disk with a continuous range in temperatures, each must contain a wide range in the actual 
temperatures. 

It is interesting to note that the tempera tures for what we call warm dust, T^;, are 



rather close to those of the "hot component" of ISridharan et al.l (|2002| ) based on IRAS data 
(Tm, of 170 K and 150 K for 18151 and 20343, respectively). It suggests that our IRTF 
measurements are measuring much of the same dust as the shorter IRAS bands measured. 
In fact, our N band F,j values for 18151 and 20343 IRS 1 + IRS 3 are each slightly more than 
50% of the IRAS PSC values. Our 24.8 /im for 18151 virtually equals the IRAS PSC 25 
fim value, and the sum of the 20343 24.8 /im F^, values is 2/3 of the IRAS 25 /xm value. 

Another interesting aspect is the extremely small amount of gas and dust in the hot 
components. This hot material is unlikely to be a major part of accretion disks that might be 
expected for HMPOs, because the accretion disks are likely to contain about the same mass 
or more as the stars, and to have characteristic temperatures from ten to several hundred K 



(jCesaroni et al.l 120071 ). It could be in a hot inner rim or surface layer of a photoevaporating 



^Of the models, x™ was largest for 20343 IRS 3. Its components, Xm,^ are 0.0005, 0.0008, 0.0035, 0.0093 
and 2.7 x 10~^ for i = 1 to 5, respectively. 
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disk (iHoUenbach et al.lll994j ). The assumed diameters of the mid-ir components suggest that 
the emission comes from material that might be described as being in an outflow cavity wall. 
In fact, it may be coming from the intersection of a flared ac cretion disk with th e surface of 
the outflow cavity, as has been suggested for other sources by iDe Buizer I (1200 7bl ). Although 
the assumed diameters of the hot components are somewhat arbitrary and affect Nn^h, they 
do not affect the mass estimates since these components are optically thin. The ranges about 
the central value of Nh^^ and are about ±30%, as expected for optically thin sources. 
(Lower limits to Nh^^ based on the 3a upper limits on source diameters shown in Table 2 can 
be calculated from the data in Table 3 since Nh is inversely proportional to Q for optically 
thin sources.) Accurate estimates of Nh^h and number density nn^h of the hot components 
of 18151 and 20343 IRS 1 will require higher resolution observations on a larger telescope. 

The optical depths in emission of the warm components, r^, are high in the cases of 
18151 and 20343 IRS 3, 4.2 and 0.9, respectively. This effect is a surprise because the 
assumed values of the projected diameters are not particularly small compared to ob- 
served and expected diarneters for structures like candidate accretio n disks near HMPOs, 
jPe Buizer fc Minierl[2005l : Ishepherd et al.l[200ll : ICesaroni et allboOTl ). While 18151 was un- 
resolved and its value for A^^,^ is essentially a lower limit, IRS 3 was resolved, so its value for 
Nh,w is a firmer estimate. For these two HMPO candidates, indicated Nh,w values are about 
as large as the column density for extinction, NH,a- Even with these high values, however, 
the masses are less than 1 Mq and much less than the mass expected for an HMPO. This 
situation again indicates that the mid-ir emitting dust is not tracing the bulk of the mass 
expected to be in accretion disks. 

Our mass estimates have led to the conclusion that in these three HMPO candidates, 
mid-ir emission is not indicating massive accretion disks. Either only a small fraction of 
their mass is emitting in the mid-ir, or the disks have been disrupted already. The emission 
may well come from d ust in and around the walls of o utflow caviti e s as has been suggested 
for other HMPOs by IPe Buizer fc Minierl koQ^ and IPe Buizer I J2006I . l2007ah and as is 
in dicated in the two dina ensional radiative transfer models for a low mass class I protostar 
of Iwhitnev et all J2003ah . 



Extinction optical depths, r^, are not extremely high, and cover a limited range. There 
is a selection effect: if they were larger, the HMPO candidates would not have been in- 
cluded in the original survey for lack of 12 /im IRAS detection, or we might not have 
det ected the HMPO candidates in the mid-ir on the IRTF. (Five of 23 fields chosen from 
the ISridharan et al.l (120021 ) survey resulted in non-detections in N band on the IRTF.) The 
appearance of the 20343 IRS 3 grism spectrum suggested that there might not be any ex- 
tinction, but fitting the spectrum required all of the components. In some cases, the range in 
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Ta from lower to upper value is small because the extinction is exponential and the sources 
are optically thick in extinction. The values of N^^a and Ay are also not extreme since 
they are directly proportional to r^. Unlike the emitting components, the extinction and 
parameters derived from it {Nfj^a and Av^a) are not affected by the assumed source size, and 
hence are more firmly defined values. 

H igh mass stars are e xpected to form at the centers of cluster-forming molecular cloud 
cores ( jBeuther et al.ll2006l ). In our small sample of three objects, only 18151 appears to be 
centered on a 1.2 mm continuum core in the plane of the sky. Comparison of the column 
density of H nucleons per cm^, Nn.a, and the visual extinction in magnitudes, Av,a, from the 
mid-ir m odel to the values for c olumn density of H2, Ngas, and Ay from 1.2 mm observations 



given by iBeuther et al.l (l2002al ) would indicate if the candidate HMPO is indeed centered 
along the line of sight within the larger 1.2 mm dust core if a consistent dust model and units 
were used at both wavelengths. If the HMPO candidate were at the core center, interior 
to the bulk of the 1.2 mm dust core's mass, we would expect the mid-ir based column 
density for the cold component to be due to near side of the dust core, and have one-half 
the column density based on the 1.2 mm observation. F or a consistent comparison, we 



have calculated from equat ion (5) using Drainel ( 
for (11^0^ at 1.2 mm given by 



nn2 

4 ^FWHM 



The 



2003al ) dust, the peak flux of 673 mJy 
Beuther et al.l fl2002al ) in their Table 2, T^ust = T^d = 47K 
given by ISridharan et al.l (|2002[) in their Table 1, and the solid angle of (11")^ in place of 
Draind ( l2003al ) Ry = 5.5 dust model has i^^e(1.2mm) = 0.2388 cm^ g~^ and 



Ay = 7.29 X where Nh is the column density of H nucleons. For 18151, the 1.2 

mm results are Nh = 5.6 x 10^^ cm~^ {hgNn = 23.75) and Ay = 410. One half of the 
column density would give logNn = 23.45 and Ay = 205. These values are a factor of 2.8 
larger than the mid-ir based values in Table 3, logiVj^^a = 23.00, and Ay^a = 72.1. It appears 
that the mid-ir source is not at the center of the 1.2 mm emitting core, but near its front 
side. However, there are significant uncertainties in the estimates. We have ignored the 
column densities of hot and warm dust in emission in the mid-infrared, but their size ~ 1" 
on the sky should not have contributed significantly to the 1.2 mm emission detected in an 
11" beam. The largest uncertainty lies in the value of Ke{1.2mm) = 0.2388 cm^ g~^. For 
far-ir through mm wavelengths, the emission coefficient is often approximated as i^' oc A^^. 
Drainel (l2003al ) Ry = 5.5 dust has P = 1.8 b e tween 250 fim and 1.2 mm. Based on a 
sample of 69 HMPO candidates, Williams et al. (2004) derive a mean /? of 0.^, that would 
increase ii'e(l-2mm) by 4.1, and decrease the column density to give logNn = 22.84 and 
Ay^a = 50.1, in much better agreement with the mid-ir based results. In fact, one could turn 



^For 18151. IWilliams et al.l ^OOi) found P = 0.5 for the limited range of 450-850 /im. However, /3 ~ 1 is 
much more commonly cited for HMPO candidates. 
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the problem around. For cases where there is strong evidence that an HMPO is centered in a 
core, the ratio of column densities determined by the extinction at 9.7 fim and 1.2 mm dust 
emission determines the overall slope of the extinction curve between them, if one assumes a 
wave length in th e far-ir at which Kg[X) flattens. Using 250 /xm as the reference wavelength 
(e.g. iHildebrandl (119831 )). the mid-ir derived column density Nh = 1.0 x 10^^ cm~^, and 1/2 
of the 1.2 mm flux density to account only for the emission in front of the HMPO, we flnd 
/3 = 1.2. 

There are other possible reasons for discrepancies. For cores that do not flU the 1.2 
mm radio beam, one would expect the difference in telescope beams to affect the ratio of 
column densities, with the larger beam 1.2 mm Nh being less than the smaller infrared beam 
Nh- Conversely, an unresolved clumpy structure in the 1.2 mm emitting core could result 
in higher 1.2 mm beam-averaged value than observed in a small diameter infrared beam 
that happens to pass between clumps giving a low extinction lin e of sight to the HMP O. 
Simil arly, viewing the source along an outflow cavity as modeled bv lwhitnev et all J2003al lbl. 
2004 ) would result in lower ir-based Nh and Ay- However, the outflows observed for these 
objects appear to lie in or near the plane of the sky (see §4). 

The HMPO c andidates 20343 I RS 1 and IRS 3 are not projected against peaks in the 
i.Z mm emission JPalau et al.lboOTh Nevertheless we can e stimate Nh and Ay for them 
from the 1.2 mm map JSeuther et al.ll2002al : [Palau et al.lboOTh that shows ~ 200 mJy/(ll")^ 



for both. The ratios of mm-based Nnf^ to ir-based Nh are 1.5 and 5.9 for IRS 1 and 3, 
respectively, in comparison to 2.8 for 18151. These are consistent with IRS 1 being deeply 
embedded but not the apparently more evolved IRS 3. 

Observations that cover both mid- and far-ir wa velengths are commonly used to estimate 
the total luminosity of an embedded source (e.g. ISridharan et al.l (120021 )). For a case in 
which an HMPO is completely embedded so that the observed infrared emission comes from 
a circumstellar dust envelope that absorbs and reradiates all of the HMPO's emission, the 
total luminosty is accurately measured by the integrated SED. Far-ir observations like IRAS 
data cover the peak of the SED, but the large beams can include sources in addition to 
that of the HMPO of interest, especially since high mass stars usually form in clusters. The 
higher resolution mid-ir observations can be limited to a single HMPO, but lack the far-ir 
contribution needed for total luminosity. Consequently, luminosity based solely on mid-ir 
fluxes is usually quoted as a lower limit to total luminosity. Our simple models can be used 
for lower limit luminosity estimates that account for extinction in the mid-ir and include 
some flux outside the mid-ir. In Table 3, we show the integrated model SED fluxes that 
would be observed through the extinction (the red curves in Figures 6-9) as We take 
these model based values as lower limits to the HMPO candidates' luminosities that are 



- 19 - 



somewhat improved over simply adding our observed mid-ir flux densties. 



In Table 3, we also show the summed modeled luminosity emitted by the hot and warm 
clouds together without the extinction factor applied (the orange curves) as If the hot 
and warm dust components had absorbed all of an HMPO's power and emitted (re-radiated) 
it as modeled, the sum of their SEDs would give the total luminosity. Presumably, the 
cold extinction layer re-radiates in the far-ir the energy it absorbed from the hot and warm 
components, so an observed full SED extending into the far-ir would have the same luminosity 
as An outflow cavity should not have much effect on the the luminosity estimate. The 
mid-ir portion of a source's SED is affected by the line of sight to t he cavity, but the far-i r 
portion that dominates the luminosity is not as strongly affected ( jWhitney et al.l l2003bl ) . 
Outflow cavities for our sources would have little effect because they appear to be close to 
the plane of the sk y (see §4). Our va lue should be close to the luminosity based on 



IRAS presented by ISridharan et al.l (120021 ). if we have observed the dominant HMPO and 
there is little luminosity from the other YSOs. For 18 151, our estimate is 22 400 Lq in close 
agreement with the IRAS-based value of 20000 L© JSridharan et al.ll2002h . The HMPO 
candidates 20343 IRS 1 and IRS 3 lie in a single IRAS beam. For 20343 IRSl and IRS3 
combined, our value of is 2200 L(7), about two-thirds of 3160 Lq from IRAS. The star in 



20343 IRS 3 is apparently more evolved than the HMPO candidate in 18151 (see §4). As a 
consequence, it has a lower column density envelope (iV^) than 18151, so its envelope does 
not absorb the full stellar emission. 

Overall, the luminosity estimates agree well with IRAS. Earlier we noted remark- 
ably close agreement between our flux densites and those of the IRAS PSC, and between 
our and the IRAS-based T^d of ISridharan et al.l (120021 ). Consequently the agreement of 
luminosities is not surprising. 



3.4. Analysis of High Resolution Spectra 

The infrared flne-structure lines are coUisionally excited forbidden lines. As such, they 
have emissivities proportional to the product of the density of the emittin g ion, the electron 



densi ty, and the collisional excitation cross section, NiNeQiu, for Ne « Nc (lOsterbrock fc Ferland 



20061). Nc, the electron density at which the collisional deexcitation rate equals the radiative 
deexcitation rate, is 5 x 10^ cm~^ for [Ne II] and 2 x 10^ cm~^ for [S IV] . The collisional 
excitation rates are proportional to 2"-i/2g--E„/fcT^ about T~°'^ for these lines. Since the 
radio free-free emissivity is proportional to NeNpT~^'^^ , the ratio of the flne-structure line 
fluxes to the free-free continuum flux is proportional to the ionic abundance relative to that 
of ionized hydrogen, with a weak dependence on electron temperature and density. The ionic 
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abundances relative to the total atomic abundances depend most strongly on the spectral 
type, or effective temperature, of the ionizing star, with weaker dependences on the stellar 
luminosity and the electron density. Consequently, the ratio of the fine-structure line fluxes 
to the free-free flux and to each other can be used to determine the stellar spectral type. 

The most convenient way to determine the stellar parameters that are c onsistent with 



the m easured fluxes is with a nebular modeling program, such as cloudy (iFerland et al. 
19981 ). We used cloudy (version 07.02.00) to calculate several dust-fre^ models with Stella, 
parameters similar to those of a B2 star, T^ff ~ 20,000 K and N^yc ~ 4 x 10^*^ photons s~^. 
An acceptable fit to the observed free-free and [Ne II] fiuxes from 20343 IRS 3 was found for 
Teff = 19,000 - 22,000 K and Niyc = 3 - 5 x 10^^ photons s-\ assuming A^'e « 5 x 10^ 
cm~^. The [S IV] fiux was predicted to be much less than that of [Nell], consistent with our 
failure to detect [SIV]. 

The extinction-corrected luminosity of IRS 1 is 1380 Lq and neither radio free-free nor 
[Nell] emission was detected from it, with limits of about 1/3 and 1/15 of the free-free and 
[Nell] fiuxes from IRS 3. The luminosity and lack of observable ionized gas are consistent 
with a spectral type of B3 or later , or Tp ff < 18 ,000 K, as the ionizing fiux of a B3 star is 



only about 1/10 that of a B2 star (|Panagialll973l ) 



4. Discussion of Individual Sources 



4.1. IRAS 18151-1208 



This HMPO candidate was chosen for study because of its deep silicate absorption in 



the mid -ir (Figure 2), the location of its mid- ir peak at a 1.2 mm peak in thelBeuther et al. 



(2002a) survey, its low l evel of 3.6 cm emissio n ( Sridharan et al. 2002 : Carral et al. Il999[ )^ 



its 

large-scale CO outfio w (IBeuther et al.ll2002br) . its detection in the K band, and the presence 
of H2 molecular jets (iDavis et al. I l2004[ ). At mm and sub-mm, there is a relatively simple 
peak 13^.''2 east and 4( ^ 9 south of the IRAS position with some extension to the southwest 
(IBeuther et al.ll2002al : IWilliams et al.ll2004l ). There is a strong second 1.2 mm peak about 
10 0^^ west and 33^^ south form the IRAS position, and there are two weaker peaks > 80" from 
it (IBeuther et al.ll2002al ). The main sub-mm peak is within ~ 2" from the MSX source that 



^ The presence of dust in an HII region affects both the free-free and [Ne II] fluxes in the same way. It 
could cause N^yc for the star to be underestimated. Evidence that N^yc were underestimated would come 
from an ir-based luminosity larger than expected from N^yc and stellar models. That does not appear to be 
the case for 20343 IRS 3 (see §4,2). 
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we identify with our mid-ir source (jWilliams et al.ll2004j ). The IRAS and m id-ir luniinosit y 
estimates agree we ll at ~ 20000 Lq (§3.3), close to that of a BO ZAMS star (|Panagialll973l ). 
Davis et al.l (120041 ) characterize the K band emission as showing a dense cluster. They 
conclude that IRS 1, the source at the K ba nd peak, is a pre- UC HII HMPO. It is the only 
deeply embedded near-ir object in the field (IDavis et al.ll2004r). a nd it is coincident with the 
mm peak and with an 0.5 mJy 3.6 cm source (jCarral et al.lll999l ). 



For an optically thin HII region, the 3.6 cm flux density of 0.5 mJy would require 
log A^i=44.67 ionizing photons s~^ at 3 kpc, using the common rel a tionshi p between flux 
density and number of ioniz ing photons as given by iDe Buizer et al.l (l2005al ). This number 
is that of a B2 ZAMS star (jPanagialll973l ). that would have L = 2900 compared to the 
ir-based luminosity of ~ 20000 Lq. The order of magnitude difference in luminosity could 
be explained by a significant contribution from lower luminosity members of the cluster, 
but this seems unlikely since the 24.8 fim diameter has a 3a upper limit of 3670 AU (Table 
2). Apparently, an HMPO at IRS 1 at the pre-UC HII stage is creating the ionization as 
a hypercompact HII region ( HC HII) or as a j et. (In either case the equation used for 
Nl does not apply.) In fact, iDavis et al.l ( l2004l ) found two lines of clumps of shocked H2 
indicating the presence of two jets, one of them centered on IRS 1, and detected Br7 at IRS 
1. Assuming Ay = 30, they estimated an accretion luminosity of ~ 120 Lq and inferred a 
total L ~ 1000 Lq by extending results from low luminosity YSOs. They argued that these 
are lower luminosity limits because their Ay was calculated for the extended H2 flow rather 
than for the source center where we have foun d Ay = 72. Our value suggests increasing 



the extinction correction for Br^ by 83 using iDraind (l2003al ) dust, resulting in an inferred 
luminosity well in excess of 20000 Lq. 



The large scale CO outflow red and blue shifted lobes are c omplex (iBeuther et al. 



2002bl ). The two jets of H2 clumps are nearly at right angles, and iDavis et al.l (120041 ) in- 
terpret the complex CO outflow as co nsistent with them. Both of the H2 jets are close to 
the plane of the sky (IDavis et al.ll2004l ). so it is unlikely that our line of sight is along a low 
extinction outflow cavity. The overall CO outflow distribution and the more powerful H2 jet 
both have their centers near IRS 1, while the s econd jet's c enter appears to be about 10" 
SW of IRS 1, well outside the mid-ir emission. iDavis et al.l (120041 ) estimate the luminosity 
of the jet centered on IRS 1 to be 0.7 Lq, and the second jet to have on ly ^ 0.05 Lp. T he 



H2 jet luminosity of 0.7 Lq is much higher than those in low mass YSOs (IDavis et al.ll2004l ). 
A larger extinction correction for the H2 luminosity, as suggested above, will increase it 
significantly. A large outflow luminosity would suggest that it is still strongly accreting. 

This source has a Class II CH3OH m aser at its mm pe ak, but neither an H2O maser, 
(IBeuther et al.ll2002d ). nor an OH maser (lEdris et al.ll2007l ). (An H2O maser is associated 
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with the second 1.2 mm source in the 18151 field.) iBeuther et al.l (l2002d ) summarize the 
conditions for models of radiative pumping of the CH3OH masers as T ~ 150 K, methanol 
column density Nm > 2 x 10^^ cm~^, and < 10® cm~^. These are are quite close to 
the parameters of the warm dust in our mid-ir emission model that has = 136K, and 
Nh = 8.5 X 10^^ cm~^ with a diameter of 2050 AU. If the source has a line of sight dimension 
equal to its diameter, the density is uh = 2.9 x 10^ cm~^. A recent set of models for CH3OH 



maser cites fractional methanol abundance Xm = 1 ^'^ (jCragg et al.l 20021 ) that would give 
methanol column density Nm = 1-7 x 10^'' cm~^. ICragg et al.l (120021 ) cite Nm/^V where 



1 km s ^ as typical. With it, our 



]^gi2.23 3 g It s couditious appear 



is the linewidth as a key parameter. They give 
model warm dust component would have Nm/^V 
to be comfortably within the conditions for strong 6.7 GHz emission (jCragg et al.l (|2002| ) 
Figure 1). Thus our mid-ir based model density and temperature are consistent with the 
observed methanol maser emission, and with a pre-UC HII stage for 18151. 



lO^^-^cm-^ 



For OH masers, the use of Xqh = 10~^ fICragg et al.ll2002r ) gives Nqh/^V 
s. This value and uh = 2.9 x 10^ cm~^ are in a region of parameter space where OH masers 
are unlikely for = 150 K (ICragg et al.l (l2002l ) Figure 1), consistent with observations. 



There have been tentative suggestions t hat masers might appear in an overlap ping 
sequence of CH3OH, H2O, and OH in HMPOs JSeuther et al. l2002cl : ICragg et al.ll2002h . but 
the evidence does not seem conclusive ( jPe Buizer et al.ll2005al ). If the sequence were correct 
it would suggest that 18151 is an early stage HMPO. Overall, our mid-ir observations and 
models strongly support the proposition that 18151 IRS 1 is an pre-UC HII HMPO, with a 
luminosity suggesting type BO. 



4.2. IRAS 20343+4129 IRS 1 and IRS 3 



IRAS 20343+4129 is one of the brighter sources in the mid-ir in the ISridharan et al. 



(I2OO2I ) list of HMP O candidates, an d its IRAS LRS spectrum shows a clear silicate absorption 
feature at 9.7 fim (IVolk et al.l 199ll). At a relatively close distance of 1.4 kpc, its IRAS based 
luminosity is only 3200 Lq (ISridharan et al.ll2002l). It has weak 3.6 cm c o ntinuum emis 



sion, bu t not any maser emission (ICarral et al.l Il999l : ISridharan et al.l |2002| ) . iBeuther et al. 



(l2002bl ) found it to have two massive molecular outfiows. The stronger massive outfiow is 
close to t he near- and mid-ir sources, but its outfiow luminosity is among the weakest they 
observed. iKumar et al.l (|2002| ) found three K band continuum YSOs with compact circum- 
stellar H2 emission. Two of them, IRS 1 and IRS 3, have the mid-ir counterparts in our 
observations (Figure 1). These two sources are oriented on an approximate NS line between 
two 1.2 mm peaks that fall on an EW line. There appears to be a partial fan of H2 emission 
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sion to be extended NE and W (see lPalau et al.l (120071 ) Fig. 1). The large scale CO outflow 



surrounding IRS 3 whose apex would be to t he north near IRS 1. ICarral et al.l (119991 ) had 
found 3.6 cm continuum at IRS 3; t he data of ISridharan et ah (2002) show the 3.6 cm emis- 



has its axis oriented NS, centered between IRS 1 and IRS 3, with red and blue lobes that 
are extended EW. 



Palau et al.l (120071 ) have observed this fleld on the SMA. They found a weak 1.3 mm 



dust peak and a C0(2-l) peak at IRS 1, with a compact EW bipolar CO outflow centered 
there. They argued that IRS 1 is ap parently a low/interrn ediate mass YSO. They also 
suggested that the redshifted CO lobe of Beuther et al. J2002bh covering IRS 1 has a different 
underlying spatial scale that the blueshifted CO lobe covering IRS 3, so that the two may 
not be directly related. Our mid-ir spectrum and 24.8 /im photometry for IRS 1 suggest an 
embedded YSO. Our model gives Ay =46 and a lumino sity Lp = 1400 Lq, consistent with a 



B3 ZAMS star that is too cool to create an HII region (jPanagialll973l ) that we could detect. 
The lack of an HII region is conflrmed by a lack of [Ne II] emission in addition to a lack of 
detected 3.6 cm emission. This object is most likely an intermediate mass YSO that is a 
significant contributer to the IRAS-based luminosity of IRAS 20343. 

The situation at IRS 3 seems ambiguous. CS (2-1) obse rvations of 20343 wit h a 27" 
beam gave a linewidth of 2.6 km s~^ centered on 11.4 km (iBeuther et al.ll2002a] ). With 
the SMA, iPalau et al.l (120071 ) detected small amounts of 1.3 mm dust emission east and 
northwest of IRS 3, and low velocity CO (2-1) emission north, east, and west of IRS 3 with 
a 3.3 km s~^ width. The large-sc ale outflow red and bl ue shifts are separated by about 5.5 
km s~^ centered on 11.5 km s~^ (IBeuther et al.ll2002bl ). but the red lobe does not appear 
to extend to IRS 3. In contrast, the [Ne II] observed FWHM is 8 km s~^, indicating an 
intrinsic width of 4-6 km s~^, depending on the line shape, wider than the molecular lines. 
It is centered on 15.7 ± 1 km s~^, close to the velocity of the redshifted low velocity (13-15 
km s~^) molecular outflow that is associa ted with IRS 1 on ly, rather than the blueshifted 
velocity ranges of either IRS 1 and IRS 3 (jPalau et al. 2007). In our analysis, it is assumed 
that the centers of the 3.6 cm emission and the [Ne II] emission are coincident. 



Palau et al.l (120071 ) suggested that the large scale massive blue shifted CO lobe is asso- 
ciated with the fan of H2 that is traced by lines of clumpy emission ~ 10" east and west of 



IRS 3 (IKumar et al.ll2002l ). They argued that the emission around IRS 3 could be powered 
by a either a B2 star with an UC HII region or a lower mass YSO with an ionized outflow. 
In either case, a cavity would have been created whose walls emit the 1.2 mm continuum 
emission condensations and the low velocity C0(2-l) se en by the SMA to the east and west, 
and the fan of 2.12 //m H2 emission (jPalau et al.l 120071 ). 



We consider the ionized outflow model flrst. iPalau et al.l (120071 ) showed that a stellar 
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wind assumed to have a 200 km s~^ velocity (jPalaij 120071 : iBeltran et al.ll200ll ) could create 
the cavity. There are NE and W extensions of the 3.6 cm emission shown in Fig. 1 of 
Palau et al.l (120071 ) that would be consistent with jets directed toward the cavity walls. The 
24.8 /im image (Fig. lb) is extended parallel to the 3.6 cm extended emission to the NE, as 
would be expected if it traced a cavity wall, and there is very faint N band emission that 
appears to be associated with the H2 emission shown in Fig. la. In contrast, the GEMINI- 
TEXES 12.8 /im continuum map (Fig. 4) is extended NS rather than NE or W that would 
be expected if it traced the inner cavity wall around the 3.6 cm emission. Compared to the 
NE-W 3.6 cm emission, the NS 12.8 fim continuum emission appears to be a possible disk 
seen edge-on, but the [Ne II] emission is at the north end rather than at the center (Fig. 5). 
If there were a 200 km s~^ wind, we might expect to see a much wider [Ne II] line than the 
observed 8 km s~^ either as part of the ionized wind or due to shock excitation. The 15.7 
km s~^ Vlsr of the [Ne II] is puzzling in this context in relation to the molecular cloud's 
velocity of 11.5 km s~^ and lack of redshifted CO emission at IRS 3. Finally, although the 
wind driven cavity hypothesis suggests the possibility of an intermediate mass protostar, th e 
lack of reddening of IRS 3 in the near-ir argues for a more evolved object (jPalau et al.ll2007l ). 
and the lack of deep silicate absorption in the mid-ir in our spectrum does as well. 



As a likely alternative to a wind driven cavity iPalau et al.l (120071 ) suggested that cavity 
could have been cleared by radiation pressure of a B2 star. The [Ne II] emission is consistent 
with the strength of the 3.6 cm emission of 1.8 mJy for the existence of an UC HII powered 
by a B2 star. The observed shift in Vlsr of about 4 km s~^ relative to the molecular cloud 
is consistent with observed dispersion of stellar velocities in clusters, and models of UC 
HII regions within molecular cloud s that sornetime s require relative velocities between the 
exciting star and its parent cloud (jZhu et al.ll2005l ). The width of the [Ne II] line of 8 km 
sec~^ suggests small velocities of expansion, stellar wind driven motion, or turbulence for an 
UC HII regionQ While the TEXES [Ne II] map does not show the NE-W extension of the 
3.6 mm continuum, the cloudy model indicates that the H"*" zone should be more extended 
than the Ne"*" zone for a B2 star. Our MIRSI mid-ir spectrum suggests an evolved envelope 
containing almost no hot dust. The extended diffuse 8-13 fim N band emission with its 
marginally resolved peak, and the marginally resolved 24.8 /im peak with relatively large 
flux are consistent with partial destruction of the inner envelope that surrounded the star, 
even though the CS, CO and [Ne II] linewidths do not suggest current high velocity or high 
luminosity outflows. The different images in the N band, the TEXES 12.8 /im continuum 



^ The Doppler thermal width for Ne in an 8000 K HII region is 4.3 km s ^. Usually H recombination 
linewidths are quoted, for which the Dopple r thermal width is 19 km s~^. UC HII regions typically show 



recombination linewidths of 30-40 km sec ^ (Hoare et al 



2007; 



Zhu et al 



200. 
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(extended NS), and the 24.8 fim filter (extended NE) suggest remnant fragments of tlie 
original disk and envelope. The remnant star format ion core has a mid-ir model luminosity 
of Le = 850 Lq, well below 3500 Lq of a B2 star (jPanagial Il973l ). but it may contribute 
more than 850 Lq to the large-scale far-ir emission detected by IRAS. This scenario that a 
B2 ZAMS powers an UC HIT that creates the 3.6 cm and [Ne II] emission seems more likely 
than that of an intermediate mass protostar creating them though a high velocity ionized 
wind . 



The high resolution mid-ir observations have identified the two most luminous sources 
in IRAS 20343+4129. Together they can account for its IRAS luminosity. 



5. Summary and Conclusions 



We have presented high resolution mid-ir observations made with MIRSI on the IRTF 
and TEXES on Gemini North of three HMPO candidates taken from a partial follow 
up survey of HMPO candi dates originally studi ed at 1.2 mm and radio wavelengths by 
Sridharan et al.l (|2002| ) and iBeuther et al.l (l2002al ). They are typical for HMPO candidates 
observed in the follow up survey being compact at 1" resolution, having low resolution spec- 
tra with strong, moderate, or weak silicate absorption, and with one emitting the [Ne II] 
line. 

A simple model of hot dust in emission, warm dust in emission, and cold dust in absorp- 
tion was developed to fit our 8-13 /im low resolution spectra and our 24.8 /im photometric 
points. Even an apparently flat 8-13 fim spectrum requires an absorption component if the 
underlying emission is assumed to be due to hot or warm silicate dust. The temperatures 
ranged from ^ 400-1000 K for the hot dust, and ~ 100-200 K for the warm dust. Using 
Draind (l2003al ) i?y=5.5 model dust properties and gas-to-dust ratio, only small masses of 
gas and dust in the two emitting components are needed to fit the data. The masses are 
less than about 1/10 solar mass (often much less) even though these are high or intermedi- 
ate mass stars, and the mid-ir emission cannot be due the the bulk of the mass in massive 
accretion disks. The mid-ir is likely to be emitted by the inner walls of outflow cavities and 
perhaps partly by the surfaces of accretion disks. On the other hand, high column densities, 
10^^ — 10^^ H nucleons cm~^, are required for the cold absorptio n components. These col- 
umn densities are less than derived from 1.2 mm 11" data using iDraind (j2003al ) dust, but 
the discrepancy may be resolved if the slope of the absorption coefficient from far-ir to mm 
is flattened, as suggested by some observations. Our three component model is not meant 
to fit either near-ir or far-ir to mm ends of SEDs. Nevertheless, the dust we are modeling 
in the hot and warm components appears to absorb the bulk of an HMPO's or intermediate 
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mass YSO's photospheric emission, so that the integrated flux of the two model components 
without application of the cold dust's extinction matches the luminosity as measured in- 
cluding the far-ir by IRAS. The mid-ir measurements together with the model thus give a 
reasonable way to determine the luminosity for individual HMPOs. 

The mid-ir emission of IRAS 18151-1208 together with weak 3.6 cm emission and other 
previous observations suggest that it is an early stage pre-UC HII HMPO whose luminosity 
is that of a BO star. 

TEXES high resolution spectra that cover emission hues from ionized gas can be used 
to determine the nature of the emission (jet or HII region) and help determine the properties 
of the underlying star. In the case of IRAS 20343+4129 IRS 1, a lack of [NcII] emission, 
a well defined compact CO outfiow, a moderately strong silicate absorption feature, and a 
dust model-based luminosity of 1400 Lq imply that it is an intermediate mass YSO whose 
luminosity is that of a B3 star. For IRAS 20343+4129 IRS 3 observed [Ne II] emission and 
3.6 cm free-free emission are consistent with a cloudy model indicating that the object is a 
B2 ZAMS star. Its weak silicate absorption and small mid-ir based luminosity suggest that 
it has already disrupted much of its natal envelope. 
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Table 1. MIRSI-IRTF Image Data 



Object 


RA{2000) 


DEC (2000) 


F^(10.5 /im)l 


F^(24.8 ^lraY 


FWHM(10.5 ^lraf 


FWHM(24.8 


IRAS 18151-1208 


18:17:58.1^ 


-12:07:25.6^ ^ 


11.2 


101.4 


I'.'Ol 


I'.'Ql 


IRAS 20343+4129 IRSl 


20:36:7.6^ 6 


41:40:08.0^ ® 


8.1 


15.6 


c.'gi 


I'.'SS 


IRAS 20343+4129 IRS3 


20:36:7.35 ^ 


41:39:52.5^ 


3.7 


86.6 


1^24 


2^41 


7 Aql 


19:46:15.6 


10:36:47.7 


72.2 


13.6 


G'.'QB 


1'.'71 


PSF Average* 










0'.'99 


1'.'83 


PSF Standard Deviation* 










0"088 


0'.'124 



^Janskys. Systematic uncertainty is ±10% for N band (10.5 ^lvc^), and ±30% for 24.8 /xm 
^FWHM of Gaussian fit to image by IRAF task imexamine 
•^3.6 cm source position J Carral et al1ll999h 
"IRTF indicated N position: 18:17:57.9 -12:07:27.0 

source position ijKumar et al.|[2002h 
''IRTF indicated N position: 20:36:07.6 41:40:10.0 
■^IRTF indicated N position: 20:36:07.3 41:39:54.0 
*Based on 6 observations of standard stars in N and 5 at 24.8 ^m 
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Table 2. Source Diameters for Models 



Object 


Distance, pc ^ 


Hot Component, AU 


Warm Component, AU 


Assumption for Calculation 


IRAS 18151-1208 


3000 


1280 


2050 


Iff deviation for unresolved source^ 


IRAS 18151-1208 


3000 


2310 


3670 


3cr upper limit for unresolved source^ 


IRAS 20343+4129 IRSl 


1400 


600 


960 


Icr deviation for unresolved source 


IRAS 20343+4129 IRSl 


1400 


1080 


1720 


3(7 upper limit for unresolved source 


IRAS 20343+4129 IRS3 


1400 


1120^ 


2200 


TEXES measurement (Hot) or 










MIRSI Gaussian deconvolution (Warm)'^ 


IRAS 20343+4129 IRS3 


1400 


1590 


2930 


3(7 upper limit for resolved source 



ISridharan et al.1 ||2002|') 

■^See text for details 

''Circular equivalent to TEXES-Gemini ellipse 



Table 3. Source Model Parameters 



Object 




Th 




\ogNH,h 


1 My, 






Tw \ogNH,w 




Ta 














AU 


K 






AU 


K 
















18151 


1280. 


455. 


0.0545 


21.05 


-3.45 


2050. 


136. 


4.20 


22.93 


-1.15 


4.90 


23.00 


72.1 


3.63 


4.35 


0.0045 


18151 


1280. 


1000. 


0.0053 


20.03 


-4.46 


2050. 


171. 


1.54 


22.49 


-1.60 


4.85 


23.00 


71.3 


3.70 


4.71 


0.0046 


20343 IRSl 


600. 


606. 


0.0055 


20.05 


-5.10 


960. 


182. 


0.12 


21.38 


-3.36 


3.09 


22.80 


45.5 


2.53 


3.14 


0.0005 


20343 IRS3 


1120. 


420. 


0.0003 


18.81 


-5.80 


2200. 


110. 


0.85 


22.24 


-1.78 


0.79 


22.21 


11.6 


2.85 


2.93 


0.0141 


Upper Values 


18151 


1280. 


1000. 


0.0691^ 


21.15 


-3.34 


2050. 


146. 


7.90'' 


23.21 


-0.88 


5.33 


23.04 


78.3 


3.73 


4.71 


NA8 


20343 IRSl 


600. 


666. 


0.0072 


20.17 


-4.98 


960. 


201. 


0.17 


21.53 


-3.21 


3.47 


22.85 


51.0 


2.58 


3.27 


NA 


20343 IRS3 


1220. 


449. 


0.0004 


18.92 


-5.69 


2200. 


112. 


1.09 


22.35 


-1.67 


1.14 


22.37 


16.8 


2.95 


3.02 


NA 


Lower Values 


18151 


1280. 


422. 


0.0036^ 


19.87 


-4.62 


2050. 


123. 


O.94IO 


22.65 


-1.44 


4.66 


22.98 


68.5 


3.52 


4.26 


NA» 


20343 IRSl 


600. 


533. 


0.0037 


19.88 


-5.27 


960. 


159. 


0.74 


21.17 


-3.57 


2.80 


22.75 


41.3 


2.40 


3.01 


NA 


20343 IRS3 


1120. 


380. 


0.0002 


18.62 


-5.99 


2200. 


105. 


0.58 


22.07 


-1.95 


0.48 


21.99 


7.1 


2.70 


2.79 


NA 



^ Assumed diameter of hot component. 

2 Assumed diameter of warm component. 

3 Optical depth of cold extinction component. 

^ Luminosity "observed" due to hot and warm dust components with extinction of cold component applied. 
^ Luminosity emitted by hot and warm dust components without application of extinction of cold component. 
6 Max Th for = 455K. 

Max Tw for Tyj = 136K. The strongest effect is at 24.8 /im. 
® Values are derived from multiple models for which the parameter's upper value shown gives Xm i — ^-^^ wavelength range most sensitive to the parameter. 

8 Min Th for Th = lOOOK. 
10 Min Th for = 171K. 
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Fig. 1.— IRAS 20343+4129 in N band centered at 10.5 fim (a) and at 24.8 /xm (b). IRS 1 is 
in the NE, and IRS 3 in the SW. The grayscale is adjusted to show extended emission. Eight 
equally spaced contours indicate the relative peak strength and extension for the central 
regions of the sources. Flux densities and FWHM are given in Table 1. The true source 
positions are most likely to be at the K band positions of IR S 1: 20''36'"7.6^ 41°40'8'.'0 (J2000) 



IRS3: 20''36'^7.3^ 41°39'52'.'5 fJ2000l flK^ar et al.l[2002[). an d the coordinates shown have 



been adjusted to fit them. The 3.6 cm source of ICarral et al.l (119991 ) is within an arcsecond 
of IRS 3. 
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Fig. 2. — Grism spectra of the three sources. 
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Fig. 3.— TEXES spectra of IRAS 20343+4129 IRS 1 and IRS 3 for [Ne II] plotted in 
velocity with respect to the local standard of rest. The [Ne II] line in IRS 3 has its center at 
Vlsr = 15.7 ± 1 km s"^ (A = 12.814 pimoru^ 780.383 cm-^ after correcting for the Earth's 
motion relative to the LSR). The gap in each spectrum at Vlsr ~ —40 km s~^ is between 
grating orders. 
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Fig. 4.— TEXES scan-mode 12.8 /xm continuum map of IRAS 20343+4129 IRS 3. The 
contours are at 1.67, 3.34, ...11.7 x 10^°Jy sr~^ The beamsize of (7/5 is shown in the lower 
right. 
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Fig. 5.— TEXES scan-mode [Ne II] (12.8 /xm) line map of IRAS 20343+4129 IRS 3. The 
contours are at 2.0, 4.0, ...10.0 x 10~^ W m~^ sr~^. The beamsize of 0'.'5 is shown in the lower 
right. 



- 38 - 



0000.0 



1) 
c 
o 



<x> 

Q 



1000.0 



100.0 



1 0.0 



1 .0 r 



RAS 18151 

Hot Dust Emission 

Warm Dust Emission 

Total Dust Emission 

Total Dust Emission with Extinotion 

Data 




10 100 
Wavelength in fvliorons 



1000 



Fig. 6. — Model individual and summed hot and warm components SEDs before extinction, 
and with extinction for IRAS 18151-1208 compared to the MIRSI-IRTF data. This model 
has a hot component temperature = 455 K. Other parameters of the model are given in 
Table 3. Emission from the hot and warm components before extinction are shown in blue 
and green, respectively, and their sum, the total dust emission inside the outer cold cloud, is 
shown in orange. The model's observed emission after extinction by the overlying cold cloud 
is shown in red, and the data are shown in black. 
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Fig. 7. — Model individual and summed hot and warm components SEDs before extinction, 
and with extinction for IRAS 20343+4129 IRS 1 compared to the MIRSI-IRTF data, as in 
Fig. 6. Parameters of the model are given in Table 3. 
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Fig. 8. — Model individual and summed hot and warm components SEDs before extinction, 
and with extinction for IRAS 20343+4129 IRS 3 compared to the MIRSI-IRTF data, as in 
Fig. 6. Parameters of the model are given in Table 3. 



